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A recent article by Wojtak et al (arXiv:1504.00718) pointed out that the local gravitational red- 
shift, despite its smallness (~ 10 —5 ), can have a noticeable (~ 1%) systematic effect on our cosmo¬ 
logical parameter measurements. The authors studied a few extended cosmological models (nonflat 
ACDM, wCDM, and ino-WaCDM) with a mock supernova data set. We repeat this calculation and 
find that the ~ 1% biases are due to strong degeneracy between cosmological parameters. When 
cosmic microwave background (CMB) data are added to break the degeneracy, the biases due to 
local gravitational redshift are negligible (< 0.1a). 


I. INTRODUCTION 


The cosmological principle, that the Universe is ho¬ 
mogeneous and isotropic on large scales, has been a 
successful theoretical basis for modern physical cosmol¬ 
ogy. In this idealized homogeneous and isotropic “back¬ 
ground Universe” we can define a cosmic redshift that 
accounts for the redshift due to background expansion 
z c = ao/a—1, where a is the scale factor inthe Friedmann- 
Lemaitre-Robertson-Walker metric, 
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and the subscript 0 denotes the quantity today. The nor¬ 
malization of a(t) is arbitrary and chosen to be ao = 1 in 
this work. Natural units c = H = 1 are used throughout 
the paper. 

We are interested in the cosmological parameters con¬ 
strained by measurements of luminosity distances. In the 
idealized background Universe, the luminosity distance 
is a function of cosmic redshift z c . Thus, the obser¬ 
vational task is to measure the luminosity distance for 
objects at different redshifts and confront the data with 
the theoretical di,{z^). 

On cosmological scales, measuring luminosity distance 
is a challenging task. Ideally, for a standard candle whose 
luminosity is known as a priori knowledge, the luminos¬ 
ity distance can be determined with the inverse-square 
law for apparent brightness. In practice, the best cosmo¬ 
logical candidate by far - the type la supernova - is not a 
perfect standard candle. The environmental dependence 
of supernova luminosity [I] and dust extinction [3] are 
accounted for by introducing some nuisance parameters 
in the supernova light-curve fitting methods (see j^j and 
the references therein for a review). Other systematics 
such as second-order corrections from gravitational lens- 
ing and directional averages have also been discussed in 
the literature mm ]• 

Redshift by its definition is a direct observable. To lin¬ 
ear order the observed redshift 2 0 b s contains three com¬ 
ponents: 0 o bs = z c + z pec + z g . The redshift due to pecu¬ 
liar velocity 2 pec adds random scatter that averages out. 
The local gravitational redshift z g accounts for the dif¬ 
ference between the gravitational potentials at the points 


of light reception and emission, which does not cancel 
if our local gravitational potential differs from the av¬ 
erage potential in the supernova environments. If we 
split the supernova samples into different redshift bins 
and consider z g for each subset of supernovae, z g can be 
redshift dependent. Moreover, in a Universe with time- 
dependent gravitational potentials, z g also contains an 
integrated Sachs-Wolfe contribution that has redshift de¬ 
pendence. Simulations show, however, that z g can be ap¬ 
proximately treated as a constant at 2 < 1 and its ampli¬ 
tude is < 4x 1CT 5 at 95% confidence level (C.L.) [7J. This 
seemingly tiny systematic error by far has been ignored 
in supernova data analysis. Ref. [7] for the first time com¬ 
puted the impact of a nonzero z g and found ~ 1% bias on 
the best-fit cosmological parameters. For modern preci¬ 
sion cosmology, 1% bias is a non-negligible effect. How¬ 
ever, for the extended models (nonflat ACDM, wCDM, 
and wo-WoCDM) studied in Ref. [7], cosmological param¬ 
eters are strongly degenerate with current supernova data 
alone. The 1% shift of best-fit parameters in the flat de¬ 
generacy direction may not be a significant effect. The 
purpose of this paper is to demonstrate that the sys¬ 
tematic errors due to z g becomes negligible when we use 
cosmic microwave background (CMB) data to break the 
degeneracy between parameters. 


II. SYSTEMATIC ERRORS WITH SUPERNOVA 
DATA ALONE 

The standard treatment of supernova likelihood is to 
absorb the random scatter ^ pec into an additional un¬ 
certainty of luminosity distance. The local gravitational 
redshift is a biased error that, in principle, cannot be 
treated in the same way. With random scatter ignored, 
a nonzero z g gives an observed redshift 

~C)bs — Z C + Zg, (2) 

and an observed luminosity distance 

dL ,obs — d^(2 c )( 1 T A>bs)/(1 T Z(f), (3) 

where d^Zc) is the theoretical luminosity distance func¬ 
tion. (See Eq. (3.3) of Ref. [7].) The extra factor 
(1 + 2 0 bs)/(l + z c ), which has been ignored in Ref. [7], 
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accounts for the fact that a local gravitational redshift 
also changes the apparent brightness via an additional 
redshift of photons and an additional time dilation. We 
find that, however, for all cases discussed in this paper 
the (1 + 2 0 bs)/(l ± z c ) correction does not lead to a no¬ 
ticeable difference. 

Following the approach in Ref. [7] , we use the redshifts 
(treated as z c ) and errors in the apparent magnitude of 
580 Union2.1 SN samples [5]. Mock data with a fidu¬ 
cial ACDM cosmology (f l m = 0.3) and a fiducial z g are 
produced with Eqs. For mock data sets with 

z g = 0, ±4 x 10 -5 , we compute the best-fit parameters, 
respectively. The 68.3%C.L. and 95.4%C.L. constraints 
on parameters are computed via Monte Carlo Markov 
chain (MCMC) simulations. 

For illustration purposes, we consider a nonflat ACDM 
model and a flat wCDM model. (See [Tj for the defini¬ 
tions of these models.) For supernova data the Hubble 
constant Hq is degenerate with the supernova absolute 
luminosity and does not affect the likelihood. The likeli¬ 
hood depends on, in the nonflat ACDM case, D& (curva¬ 
ture of the Universe) and D m (fractional energy of mat¬ 
ter), and in the roCDM case, w (dark energy equation of 
state) and D m . The dark energy fractional energy Da is 
constrained by D*, + D m + Da = 1 (Dj, = 0 in wCDM 
case). 

As shown in Fig [TJ we find that the systematic errors 
on cosmological parameters agree with Ref. [7j. However, 
the ~ 1% bias of best-fit parameters are aligned with 
the degeneracy direction where the posterior likelihood 
is flat. Compared to the ler (68.3% C.L.) constraint, the 
systematic errors are small. 


III. SYSTEMATIC ERRORS WITH 
SUPERNOVA + CMB 



FIG. 1. The systematic errors due to a local gravitational 
redshift. Upper panel: nonflat ACDM model; lower panel: 
flat u>CDM model. The contours are 68.3% C.L. and 95.4% 
C.L. constraints for the z g = 0 case. 


In this section we combine the publicly available 
Planck CMB likelihood [9( and the Union2.1 supernova 
likelihood to constrain cosmological models, with differ¬ 
ent fiducial z g assumed. The Planck likelihood includes 
CamSpec temperature power spectrum likelihood (2013 
release) and a low-f? WMAP polarization constraint. The 
supernova likelihood code is modified to include the z g 
correction. 

We run Monte Carlo Markov chains with the publicly 
available code COSMOMC [10], with the standard set¬ 
tings, namely, flat priors on six parameters fl^h 2 (baryon 
density), D c h 2 (cold dark matter density), 9 (angular ex¬ 
tension of sound horizon on the last scattering surface), 
r (CMB optical depth), lnA s and n s (logarithmic am¬ 
plitude and index of the primordial power spectrum of 
curvature fluctuations). For nonflat ACDM model and 
wCDM model we add D^ and w with a flat prior, respec¬ 
tively. 

As shown in Fig. [2] the impact of z g becomes invisible 
with the addition ofCMB data. In Table |I| and Table HT1 
we compare the mean and standard deviations of param- 


TABLE I. Nonflat ACDM; CMB + SN 



Zg =0 

z g = - 4 x 10~ 5 

Zg = 4 x 10“ 5 


0.344 ±0.042 

0.344 ± 0.043 

0.345 ± 0.042 

Da 

0.665 ±0.033 

0.665 ±0.032 

0.664 ±0.032 

Dfc 

-0.009 ±0.011 

-0.009 ±0.011 

-0.009 ±0.011 


eters and find that in all cases the biases due to a nonzero 
z g are less than 0.1 <7. 


IV. CONCLUSIONS AND DISCUSSIONS 

The findings of Ref. [7], that a local gravitational red¬ 
shift can cause ~ 1% systematic errors, raise a question 
whether or not we should include the z g correction in the 
supernova data analysis. In this paper we have demon¬ 
strated that, for current supernova data, as long as CMB 
data are included to break the parameter degeneracy, the 
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FIG. 2. The 68.3% C.L. and 95.4% C.L. constraints on nonflat ACDM cosmology (left panel) and wCDM cosmology (right 
panel). The red, blue and black lines are for z g = 4 x 10 _s , z g = 0, and z g = —4 x 10 -5 respectively. Data sets: Planck 2013 
CamSpec + WMAP low-£ polarization ± Union2.1 Supernova. 


TABLE II. Flat wCDM; CMB ± SN 

Zg = 0 _ Zg = -4 X 10~ 5 Zg = 4 X 10~ 5 

Q. m 0.293 ±0.023 0.292 ± 0.023 0.294 ±0.023 

w -1.093 ±0.087 -1.096 ± 0.086 -1.089 ±0.087 


impact of the local gravitational redshift is negligible. 
Thus, no modification of current supernova likelihood is 
on demand. 

In Sec. H] and Sec. [ill] we have used different sets of 
parameters. In particular, in the supernova-data-alone 
case we have used a flat prior on Q rn . and for CMB ± SN 
we have used flat priors on f \h 1 2 3 4 5 , tt c h 2 and 0. (I2 m is a 
derived parameter with a complicated prior.) To make a 


more direct comparison, we replace the CMB likelihood 
with a “compressed CMB likelihood” that contains only 
the distance information, that is, a Gaussian constraint 
R = VnrnH 0 d A {z*) = 1.7488 ± 0.0074 [Hj, where cIa is 
the comoving angular diameter distance and z* = 1089 
is the redshift at recombination. The advantage of using 
this compressed likelihood is that R can be written as a 
function of f l m , Qk and w. For supernova and this com¬ 
pressed CMB likelihood, we use the same parametriza- 
tion and priors as in Sec. [IT] and find that the systematic 
errors due to z g remain negligible (< O.lcr). 

Whether a systematic error can be ignored depends on 
the precision of observations. Local gravitational redshift 
may still be a potential source of information contami¬ 
nation for future supernova observations. We leave the 
forecast of forthcoming supernova data to future work. 
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